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Abstract 

The Davies-Fulhng (DF) model describes the scattering of a massless field by a non-inertial 
mirror in two dimensions. In this paper, we generalize this model in two different ways. 
First, we consider partially reflecting mirrors. We show that the Bogoliubov coeflicients 
relating inertial modes can be expressed in terms of the frequency dependent reflection 
factor which is specifled in the rest frame of the mirror and the transformation from the 
inertial modes to the modes at rest with respect to the mirror. In this perspective, the 
DF model is simply the limiting case when this factor is unity for all frequencies. In 
the second part, we introduce an alternative model which is based on self- interactions 
described by an action principle. When the coupling is constant, this model can be solved 
exactly and gives rise to a partially reflecting mirror. The usefulness of this dynamical 
model lies in the possibility of switching off the coupling between the mirror and the fleld. 
This allows to obtain regularized expressions for the fluxes in situations where they are 
singular when using the DF model. Two examples are considered. The first concerns 
the fiux induced by the disappearance of the refiection condition, a situation which bears 
some analogies with the end of the evaporation of a black hole. The second case concerns 
the fiux emitted by a uniformly accelerated mirror. 
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1 Introduction 



The Davies-Fulling (DF) model [|l| describes the scattering of a massless field by a non- 
inertial mirror in two dimensions. It has received and continues to receive attention 
0, 1^, 1^, ^, H, ^ [l^, [ll|] principally because of its simplicity and its relationship with 
Hawking radiation [Q. (One can indeed mimic the varying Doppler effect induced by the 
collapse of a star by the total reflection on a receding mirror.) Because of its simplicity, this 
model has been also used to investigate the relationships between the particle description 
of fluxes emitted by the mirror and its fleld description based on the two-point Green 
function. The motivation behind this analysis is the following. When quantizing a fleld 
in a curved space-time one looses the uniqueness of choice for the particle notion which is 
then used to deflne the vacuum and its excitations. Based on this fact, some authors have 



proposed to discard the particle point of view[|l^]. The DF model, being deflned in flat 
space time and giving rise to particle creation as in a curved space-time, provides a good 
playground for confronting the two points of view. Finally, the DF model also provides 
a good starting point for studying the role of ultra-high frequencies which arise in the 
presence of event horizons [0, |T3], |TB|, [T^. This is particularly true when considering 



uniformly accelerated mirrors p, |T8|, |T9[. Indeed, in this case one has to confront the fact 



that the instantaneous value of the energy flux identically vanishes whereas the Bogoliubov 
coefficients, mixing positive and negative frequencies, do not vanish and lead to a total 
energy which furthermore diverges. 

Quite independently of these specific difficulties, there is a fundamental reason which 
renders the analysis of these problems complicated: the DF model does not follow from 
an action principle. In fact, the refiection condition is imposed from the outset instead 
of following from interactions with the boundary. Therefore only questions concerning 
asymptotic properties of asymptotically inertial mirrors can be properly answered. To 
emphasize this point, we shall show in the first part of this article that the scattering 
in the DF model can be expressed in purely kinematic terms. It results from the Bo- 
goliubov transformation relating the usual Minkowski modes to non inertial modes which 
are eigen-modes with respect to the proper time of the mirror. The scattering of the 
latter is then trivial, as trivial as the scattering of Minkowski modes by an inertial mir- 
ror. This rephrasing of the DF model is very useful in that it allows to consider partially 
transmitting mirrors with arbitrary frequency dependent transmission coefficients. In this 
perspective, the DF model is simply the limiting case in which the refiection is total for 
all frequencies. 

In the second part of the paper we analyze an alternative model for the scattering 
along a given trajectory which is based on self-interactions described by an action prin- 
ciple. The main motivation for considering this model is that we can now switch on and 
off the coupling between the mirror and the field. Therefore, we can work with well de- 
fined asymptotic free states. The relationship between the partially transmitting mirrors 
previously considered and this model based on asymptotic states will be explicitly made. 

In the first part of this Section we work with a coupling which is constant. In this 
case, the Born series can be exactly summed and lead to a partially transmitting mirror. 
Moreover, in the large coupling constant limit, one recovers the DF model, i.e. total 
refiection. Secondly, we consider the case when the coupling is time dependent. In this 
case, we compute the fiuxes perturbatively to quadratic order in the coupling. The novelty 
arises from transients effects associated with the switching on and off. The possibility of 
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controlling these transients is crucial for regularizing the fluxes in situations where they 
are singular when using the DF model. 

To make this explicit, we consider two examples. The first one consists in computing 
the flux associated with the appearance (or disappearance) of the reflecting boundary 



condition. This problem was considered by Anderson and deWitt Moreover, as 

discussed in 0, it presents some analogies with the residual flux associated with the 
disappearance of a black hole at the end of the evaporation process. When using the DF 
model, the flux is singular and its spectral properties are ill defined. On the contrary, 
with the self-interacting model, it can be described by a well-defined process in which 
the coupling to the mirror is switched off more and more rapidly. The second application 
concerns the flux emitted by a uniformly accelerated mirror. In the DF model, the energy 
flux vanishes everywhere but on the horizons where it is not defined. With the other 
model instead, a well-defined and regular expression will be obtained. In the intermediate 
regime we recover the vanishing of the local flux. But we also find transient effects which 
become singular when the switching on and off of the coupling is performed for asymptotic 
early and late proper times. 

Finally, we wish to stress that in this paper the recoil effects shall be totally ignored 
since the trajectory of the mirror is given once for all. Nevertheless, since the interacting 
model is based on Feynman diagrams, it prepares the analysis of taking into account 
the dynamics of the mirror [pi]|, rT|. Indeed the 5'-matrices computed with or without 



back-reaction effects possess a very similar structure. 



2 The kinematic models 

In the first part of this Section, we review the basic properties of the Davies-FuUing 
model. In particular we compare the particle description of the fluxes based on Bogoliubov 
coefficients with that based on the two-point functions. In the second part, we show 
how the scattering process can be generalized so as to describe partially transmitting 
mirrors. This generalization will be performed in a matrix formalism. We have chosen 
this formalism for two reasons: first to emphasize the kinematic nature of the DF model, 
and secondly to introduce in natural terms the generalization to partial reflection. In the 
third part, we relate the Bogoliubov coefficients to the S*- matrix acting in the Fock space, 
thereby preparing the analysis of transition amplitudes performed in the next Section. 



2.1 The Davies- Fulling model 

In the Davies-Fulling model, the mirror is perfectly reflecting for all frequencies and its 
trajectory is chosen from the outset. Moreover, no width is attributed to the reflecting 
condition, i.e. it acts like a delta in space. Beside the fact that the trajectory is always 
time-like, we shall also impose that it is asymptotically inertial. In conformal terms this 
means that the trajectory starts from i~ and ends in i+, the past and future time-like 
infinities respectively The reason is that in the other cases, i.e. when the mirror 

originates and/or ends on null infinities, the calculation of the energy radiated by the 
mirror is ill-defined. (The specific problems associated with such asymptotic trajectories 
will be considered in a next article ||22|| ) . 



In this paper, we shall consider the scattering of a complex massless scalar field. The 
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reason for this choice is that it allows to consider more general scattering matrices when 
the reflection condition is not perfect. This possibility will be exploited in the next 
subsections. Since the dynamics of the mirror is trivial, the evolution of the field is 
entirely governed by the d'Alembert equation 



and the reflection condition 



{dl-dlm,z) = Q (1) 



$(t,z,K^)) =0 (2) 



along the classical trajectory of the mirror z = Zci{t). 

Since the field is massless and since we work in two dimensions, it is particularly useful 
to work in the light-like coordinates defined hj U,V = t ^ z. For instance, the general 
solution of eq. (|I]) is the sum of a function of U alone plus a function of V. Moreover, since 
the mirror is perfectly reflecting, the trajectory of the mirror completely decouples the left 
hand side configurations from the right hand side ones. Therefore, in this subsection, we 
can and shall restrict our attention to the configurations living on the left of the mirror. 

Finally, since the mirror trajectory emerges from i~ , V = —oo is a complete Cauchy 
surface. Hence, the usual modes of the d'Alembertian given by 



Q—ikU 



MU) = (3) 

/47r A; 



form a complete and orthonormal basis. (Instead, when the trajectory starts from the left 
part of J'~ , the choice of a complete and orthonormal basis should be reconsidered [p^] . ) 
We recall that the norm of the modes is determined by the Klein-Gordon scalar product 
which reads, when evaluated on J'~ : 

/ + 00 ^ 
dUipliduifik' = sgn{k) 6{k - k') . (4) 

The scattered mode corresponding to eq. (H) is determined by eq. @) to be 

^r\v) = — ^= , (5) 



where U = Uci{V) is the trajectory of the mirror in the light- like coordinates. 

The m-mode ^p^iP'{U, V) is by definition the solution of eqs. (|I|) and (0) which has eq. 
@ as initial data. It is given by 



f,-ikU ^-ikU,i(V) 



vTiu,v) = ^=- . (6) 



To analyze the frequency content of its scattered part, it should be Fourier decomposed 
on the final Cauchy surface U = +oo (the left part of J^~^)- In total analogy with what 
we have on J'~ , on J''^ the normalized modes are 



^.{V) = -^== . (7) 
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Then the scattered mode can be decomposed as 



where the coefficients a^ik, Puik are given by the overlaps 

r+oo pi^^V ^-ikUciiV) 



47r|u;|-i y47r|A;| 

p:k = {^:\^r')- (9) 

Since both the initial and the final set of modes are complete, the coefficients a^^kyPujk 
satisfy the relations 

dk {alf^a^'k - = S{uj - uj') 

duj {a^kalk' - (3ujk(3*k') = ^{k - k') 
dk {a^kf^wk' - f^ujkacjk') = 

duj {a^kPl,^ - (3lka^ik) = . (10) 



















Note that these relations are not trivially fulfilled when the trajectory of the mirror reaches 
one of the null infinities rather than the time-like ones. Notice also that the overlaps (j^) 
can be computed on any space-like surface which runs from z = — oo to some point on 
the mirror [t, Zci{t)). In this case, one should use the full expression of the in modes given 
in eq. (P) as well as that of the out modes given by 

^:'^\u,v) = ^uv)+v'r\u) . (11) 

The second term yj^^'^"* results from the backward scattering of ip^j given in eq. (^. 

When the overlaps and jS^^k are known, the classical scattering problem is solved. 
That is, it suffices to decompose the initial data in terms of the modes (|^ to obtain, 
through eq. (^), the Fourier content of its image on J'^. It should be pointed out that 
the coefficients (3i^k which mix positive and negative frequencies have a well defined role in 
this classical wave theory: they determine the (non-adiabatic |2^) increase of the Fourier 



components of the scattered wave (see e.g. eq. (11) in for their infiuence on the energy 



of the refiected wave). It should be also pointed out that one can recover an approximate 
space-time description of the scattering when considering localized wave packets rather 
than plane waves: for sufficiently high frequencies (i.e. short wave lengths compared to 
the (acceleration of the mirror)"^), the coefficients /3^fc vanish and the mean frequency 
of the refiected packet u is related to k, that of the incident one, by the Doppler effect 
UJ = k dvUci\u=u evaluated at U, the retarded time of the center of the incident packet. 
These two properties are easily obtained by evaluating the integrals in eq. (^ by the 
saddle point method. 

When a^k and Pi_^k are known, the quantum scattering problem is also solved. This 
follows from the linearity of eq. (|l|) and eq. (|^): when working in a second quantized 
framework, being linear, these equations provide the Heisenberg equations for the field 
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operator. Thus the field operator can be written both in the in and the out basis by 



[(^k V^fc + Ok Vk 



:i2) 



When imposing that it satisfies the equal time commutation relation [(^{z),dt^\z')\ = 
i6{z — z'), eq. (^) guarantees that the m-operators ak,bk satisfy the usual commutation 
relations leading to the particle interpretation. Then, the in vacuum |0j„) is defined as 
the product of the ground states of the in oscillators and its excitations are generated by 



the creation operators a^" ' , ^aT^- Moreover, by construction of the in modes on J'~ , the in 
particles correspond to the usual Minkowski particles on J'~ . Similarly, by construction 
of the out modes, all these properties apply to the out operators and to the 

out vacuum \Oout) when replacing J'~ by J''^. 

Given the orthonormal and complete character of the in and out mode basis, eqs. (j^) 
and (|l^) determine the Bogoliubov relations : 



m 
"A: 



J 

DO 



= / dk (a 
Jo 



Luk "fc" f^-^k ''') 



(13) 



Then eqs. ( pil| ) guarantee the compatibility of the particle interpretation in each basis, i.e., 
both in and out operators obey the canonical commutations relations. With the relations 
(P^), all questions concerning quantum scattering processes can be answered. For instance, 
the probability amplitude to obtain a given final state \^ fin) specified on in terms of 
out operators starting from some in state |Hj„) constructed on J~ is given by the (Fock 
space) product More intrinsic is the overlap Z^^ = {Ooutl^in) between the 

initial and final vacuum states. Indeed, it determines the probability amplitude for the 
(spontaneous) decay of the vacuum. The computation of Z is easy when the scattering 
is stationary, see e.g. @, H. In the general case however, due to the frequency mixing 
between in and out modes, the calculation of Z is less trivial. This generalization is 
presented in Appendix A. 

It should also be noted that the Bogoliubov coefficients themselves are given by the 
following matrix elements 



a,. 



{0^ 



0^n) 



-(^tk = (o,,„ir*^ar^ioj . 



(14) 



However it is not clear how to attribute a physical meaning to these equations. In partic- 
ular, the relationship with the second one and pair creation amplitude is quite obscure. 
Indeed the probability amplitude to obtain on J7+ one pair of quanta of frequencies u and 
u' in the in vacuum is given by 



1 r°° 

{Oout\a°jXf'\Orn) = dkP, 

Z Jo 



'uk^ku' 



(15) 



We shall return to these questions of interpretation in subsection 2.3. 

Instead of considering in — out matrix elements in Fock space, more attention has 
been put on expectation values of (local) operators in a given initial state. The most 
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studied object is probably the energy flux emitted by the mirror when the state of the 
field is in vacuum. The motivations for this analysis are, first, its relevance for black 
hole radiation||2[-|Tl|], secondly, that its non- vanishing value is due to spontaneous pair 
creation, a specific feature of quantum field theory (QFT), and thirdly that this value can 
be computed both from using eqs. (|T3[) and from the properties of the Green function of 
the field. 

Having at our disposal the Bogoliubov coefficients a^kyPtuk, we start with the particle 
point of view. We consider the density energy of the emitted flux. The corresponding 
hermitian operator is^ Tyy = + dv^dy^^ ■ On the left of the mirror {U > 

UciiV)), using eq. (|13D and the first line of eq. (p!OD, the expectation value of the energy 
flux is 

(Tyy) = (Oj„|Ti/y |0j„) — {Oout\Tvv\Oout) 



2 3?e 



27r 



(/"^^ - e-^(-'+-)^ (J'^dk al^f3^,k 



(16) 



It should be noted that the subtraction of the out vacuum flux follows from the prescription 
of subtracting the contribution of the Minkowski vacuum. Indeed, by construction of the 
out modes, they coincide with the usual Minkowski modes on . 

The total energy emitted to is obtained from integrating (Tyy) over V . The 
integration eliminates the second term which is due to interferences between states with 
different particle numbers. It gives 



/ + 00 
dV (Tyy) 
-oo 



-oo 

2 / duj LJ dk \Pujkf = 2 / du uj {nj) . (17) 
Jo Jo Jo 



One gets the usual relationship between the mean energy and the mean number of particles 
{nj) = J^dk IPi^kf found on J'~^ (it equals to the number of anti-particles). In this writing 
one sees that the non- vanishing character of (Hy) is due to the (3 coefficients which govern 
the vacuum decay, see eq. ( p.09| ) in App. A. 

We now reconsider the flux (Tyy) without making use of the Bogoliubov coefficients 
and with less emphasis on the notion of particle. This method is based on the two-point 
Wightman function evaluated in the in vacuum 



oo 



(o..|<i>t(f/,v)$(f/',r)|o.„) = / dk^Tiu,v)^riu',v') . (18) 

Jo 

In terms of this function, using eq. (BI), the mean flux on J'^ reads 



(Tyy) = 2]im^[{0in\dv^^dv'^\0in) - {Oout\dv^^dv'^\0out) 

= lim dvdv [\n\U,iiV') - U,i{V)\ - \n\V' - V 



^ The symmetrization is due to the fact that we deal with a complex field. Of course, in the DF 
model, particles and anti-particles equally contribute to (Tyy). This explains the overall factors of 2 in 
the next equations. We warn the reader that this equal contribution will not be necessarily found when 
considering partially transmitting mirrors. 
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1 
1 

24^ 



dV 



1/2 



d 



dV^ 



2 ( / 


'dUd 




dV 


dUci 


r\ 


dV 





-1/2N 



127r 



V 




(19) 



Again, the subtraction of the out vacuum flux follows from the prescription of subtracting 
the contribution of the Minkowski vacuum. In this second description of the flux, it is 
through this prescription that the notion of vacuum decay is re-introduced. Indeed, on 
JT""*", the above subtraction is equivalent to normal ordering with respect to out operators. 
(This is straightforwardly proven by using eq. (p!OD.) Moreover, the fact that {Tyv) 
identically vanishes for inertial trajectories, i.e. when dyUd = 0, confirms that its non- 



vanishing character is due to the non-adiabaticityP5[| of the scattering, a notion deeply 
rooted to the spontaneous creation of pair of particles. 

^From eq. (|T^) we learn that the energy flux is local in that it depends only on three 
derivatives of the trajectory UdiV) evaluated at the retarded time V (remember that we 
are on the left of the mirror). We shall see below that this locality is a consequence of 
dealing with a perfectly reflecting mirror for all frequencies. 

Notice finally that in eq. (0), the first term is positive definite and leads to the 
positive total energy {Hy)- Indeed, being a total derivative, the second term does not 
contribute to {Hy) when the trajectory is asymptotically inertial since UdiY) ~ V for 
asymptotically late and early Vs. This might not be the case for trajectories which 
enter or leave the space through the null infinities, because of the infinite Doppler effect 
encountered asymptotically. 



2.2 Partially transmitting mirrors 

In this subsection we study partially transmitting (but still recoil-less) non-inertial mir- 
rors. We shall proceed in three steps. We first show that the scattering by a non-inertial 
mirror is most simply described in terms of the wave functions which are eigen-modes 
of the proper time of the mirror. (We shall call them the proper-time modes.) When 
using these modes, the matrix relating the scattered modes to the initial ones is diagonal 
in the frequency, exactly like for the scattering of Minkowski modes by a mirror at rest. 
Secondly, we shall see that these modes are well adapted to introduce partially reflecting 
coefficients with arbitrary frequency dependent phase and amplitude. Indeed, since this 
matrix is diagonal in the proper-time frequency, unitarity constraints its elements in a 
simple manner, frequency by frequency. Thirdly, for both partial and total reflection, we 
shall see that the usual Bogoliubov coefficients, eq. (^), relating the in and out Minkowski 
modes are simply obtained from this diagonal matrix. 

To fulfill this program, we first need to construct the proper-time modes. To this 
end it is very useful to introduce new light-like coordinates u, v such that the time-like 
coordinate {u + v)/2 = r is the proper time of the mirror and the space-like one defined by 
{v — u)/2 = p is such that the trajectory reads p = po = constant. These new coordinates 
are defined by two analytic functions u{U) and v{V) where U, V are the Minkowski light- 
like coordinates. These functions are determined by the mirror trajectory Uci{V) and the 
two conditions defining r and p. Indeed, along the mirror's trajectory, the length element 
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obeys 



This gives 



= dv^ = du^ = dr^ . (20) 



dv rrr— , du 



IdvU^i and — = ^duV,i. (21) 



One verifies that the Jacobians dv/dV and du/dU define a time dependent boost since 
they satisfy {dv/dV){du/dU) = 1 for all r. The proper-time modes are then simply given 
by 



^—iXu 



47r|A| 



Vx{v) = -r=- (22) 
/47r A 



They form a complete basis on J'~ and J'^ since our condition that the trajectory emerges 
from i~ and finishes on i"*" implies that the v and u axes cover those of V and U respectively. 

In the case of total refiection, the scattering along the mirror at p = po is trivial. 
When using the conventions of the former subsection (eqs. (|^) and (0)), one has, on the 
left of the mirror, 

= -e^'^^po ^l''>^\u,v) . (23) 

The new subscript U, V indicates which side of J'~ {J^) is the asymptotic support of 
the in {out) functions. We have introduced it in order to describe partial refiection which 
requires to consider simultaneously both sides of the mirror. Using this notation, on the 
right of the mirror, one has 

= -e~2*^^V^'°"*(M,t;) . (24) 
It will be useful to express these relations by a 2 x 2 matrix Sa as 

i,out Qij j,in f j,in Qtj i,out 



= swr ( = ar = s^iar) . (25) 

At fixed A, the indices of rows and columns j are the new subscript U or V. As usual, 
repeated indices are summed over. For total refiection, one has 



— g+2iApo 



(26) 



We now consider partial refiection. When considering elastic refiection, the matrix Sa 
relating in and out modes which generalizes eq. (^) is unitary. (That is, we generalize 
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total reflection in a restricted way since we keep both the linearity and the unitarity of 
eq. (H).) Unitarity constraints the elements of Sa 



Su e''^^ -i R e'^ 



-i R' e'^ e*^" 



(27) 



to obey 

R = R! 



1 



Su + R =1 and ip' = ipu + ip^ - (f . (28) 

(For simplicity of the expressions, we have not written the argument A but all variables 
should be understood as A dependent.) Physically, R and s correspond to the reflec- 
tion and transmission coefficients, i.e. when working in the rest frame of the mirror, the 
probability for an incident quantum of frequency A to be reflected is R^. 

In what follows we impose ipu = fv = —<Pi a condition which expresses that the 
transmitted part of scattering is independent of the sign of the momentum. In anticipation 
to Section 3, we point out that this equality is automatically satisfied when considering 
parity invariant hamiltonians (see |2i| Chap. 3.4). In this case the matrix reads 



^^-"^ \ -iRe-^' ) ■ ^ ^ 

In principle the common phase e~^'^ could be re-absorbed in a redefinition of the modes. 
However, when using in and out modes conventionally defined, i.e. ip^'^"'{v) = (p\'°^^{v) = 
ipx{v) of eq. (0), the phase is univocally fixed. As we shall see in the next Section, 
this convention is automatically used when considering interactions perturbatively. This 
is also the case in the DF model. Indeed the limiting case of total reflection given in eq. 
(^) is reached for i? — >■ 1 and (f) = n /2 for all A. One also finds that the other phase 6 is 
related to the mirror location by 6' = 2poA. 

To complete our second step, we should describe particles and anti-particles simulta- 

1 i-r j-U- J j-t, • J- / tj.in V.in lU.in] iy,in]\ ■ a , 

neously. io this end, we group the m operators {a^^ ^a^ ,0^ 'j m a 4- vector 

u,m 1 J.1 J. J. / U.out V.out lU.out] TV,out]\ • ii.out o- •! i j-i • 

03^ and the out operators [a^^ , a^^ ,0^ , b^^ ') m . Similarly, we group their 
corresponding modes in the 4-vectors v^^'™ and <^^'°"*. Since we work with a charged field, 
the modes associated with might not be the complex conjugates of those associated 
with a\. (As it is the case when dealing with a charged field in an electro-magnetic field, 
see e.g. Sect. 1.3 in |g]). Explicitly, in our case, the 4 modes are {ip^,ip\,ip^* ,ip\*) 
where ip\ designate the two modes associated with anti-particles operators 6^ and h\. 
We then introduce the 4x4 matrix given by 

Saa' = 5(A - A') { g° \ . (30) 



A 



where Sa is the scattering matrix for the anti-particles. Since Saa' is block diagonal, 
unitarity constrains Sa and Sa separately. Saa' acts on the in 4-vector as follows 



= Sxx' ^x' (31) 

where continuous repeated indices are integrated from — *• 00 and discrete ones summed 
over the four components defined at fixed frequency. With these choices, the components 
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of Saa' are the Bogoliubov coefficients conventionally defined. By conventionally defined 
we mean the equations which generahze eq. (plSf ), i.e., 



j.m 
Oa' 

lJ,in t 



11 
"aA' 



i,out 
«A 



nij * i,out 

P\y ^A 



+ a 



AA' 



ii,out t 
ri.out t 



(32) 



where the Bogoliubov coefficients a, a, /3, /3 are now 2x2 matrices. By direct identification, 
one obtains 



"aa' 


/ j,in\ i,out\ 


- '-'xy ' 


"aa' 


1 —i,m\ —i,out\ 

= Wy ) = 


Qi+2j+2 * 

- '-'xy 


^% 


= Wy ) 


qij+2* 

— '^xy 




1 j,in\ —i,out *\ 

= iVy Wx ) 


Qi+2j 

— '-'xy ■ 



(33) 



When Saa' is block-diagonal in the sense of eq. (^ 
In full generality, Saa' satisfies unitarity in the following sense 



one obviously has [3^-' 



xy 



Pi' 



(^t) 



XX" '-'A"A' 



6{\ - A') 6""' 



(34) 



This equation generalizes eqs. (p^ ) to partially transmitting mirrors. 

With eqs. (p6|), (|30|) and (|3TD, we have shown that the scattering in the DF model is 
trivial when using the proper-time modes. We have done more since eqs. ( PT| ) and ( P3| ) 
apply to all partially transmitting mirrors governed by Sa given by eq. (p9|). 

The last step consists in finding the relationship between Saa' and the Bogoliubov 
coefficients between in and out Minkowski modes. This is simply achieved by introducing 
the 4x4 matrix which relates the (unscattered) Minkowski modes of frequency k = —idt 
to the (unscattered) proper time modes of frequency A = —id-r : 



(35) 



The elements of this matrix are given by 



kX 



1 {^"iV 


fk) 









fk) 














fk) 











f) 















f*) 







V 

















f*) I 

























"fcA 





PkX 


* 










aUU 
PkX 





"fcA 















V 


aVV 
PkX 





"fcA 


* ) 











(36) 



Since Bkx relates unscattered modes, it is independent of the charge of the particle, hence 
Bll = oil'i = a^A^ = Bll The same equality applies to al^ , , s^nd (3^!^. 

The important point for us is that Bkx also relates the in Minkowski modes to the in 
proper-time modes and the out Minkowski modes to the out proper-time modes. Therefore 
the linear relation between in and out Minkowski modes is given by 



LfM,OUt 



'-'ujk 



1^)1 ,in 



(37) 
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where 



Repeated indices are summed over and the inverse of B is defined by 



kX 



(38) 



(39) 



It is given by 



B-' 



kX 



1 


uu* 

"fcA 





aUU* 
PkX 





\ 





vv* 

"fcA 





aVV* 
PkX 




PkX 





"fcA 







V 





PkX 





VV 
"fcA 


/ 



(40) 



Exphcitly, using the dictionary (Rl) now apphed to Minkowski modes, the four coefficients 



^wk 



a 



uu 

uik 



5{uj — k) — i 



rl\ f^,UUrpUU^UU * , aUU *rpUU * aUU\ 
"^(«a;A^A «fcA + PloX ^A PkX ) 



—I 



^\ (^UUrpUV^VV* I oUU *rpUV * nVV\ 
(a<^A -^A "fcA + Pu^X -?A PkX ) 



UU* 
'ujk 



^\ ( ^UUrpUU oUU * I nUU *rjnUU *^UU\ 

("..A -^A PkX + Pux a^x ) 



uv* 

'u)k 



UUrpUV oVV * 



PL' * + Ptx 



uu *rpuv * „yv 



X 



a 



kX 



(41) 



Similar equations give the expressions for the remaining components of S*^^. We have 
written Sa as Sa = 1 — « Ta (and Sa = 1 — i Ta) in order to extract the trivial part of the 
diagonal elements. This trivial part leads to the delta function in the first equation. The 
usefulness of the writing is that it will be easily related to the perturbative expressions 
we shall encounter in the next Section. 

Equations (41) are the central result of this Section. They give the in — out overlaps 
of Minkowski modes in terms of the matrices Ta, Ta computed in the rest frame of the 
mirror and the overlaps between the free (unscattered) Minkowski and proper time modes. 

It is then easy to obtain the mean fiux emitted by this partially transmitting non- 
inertial mirror when the initial state of the field is the Minkowski vacuum. The same 
algebra which gave eq. (ITB) now gives 



{Tvv) = {Tvv) 



particle 



+ {Tvv) 



anti— particle 



(42) 



where 

irji \ particle 

\J-vv) 



E 

j=U,V 



,\/ Ijjijj' 



dijjdjjj 
27r 



-i{uj'-<jj)V 



d^ P^k f^^'i 



'i{uj'+Uj)V 



f 

Jo 



dk all ^ 



Pul'k 



(43) 



(Tvv) anti-particle giveu by the same expression with a,P replaced by a, (3. {Tyv) pos- 
sesses the same structure as eq. (|T6|). However, four kinds of coefficients should be 
considered since we are dealing with partial refiection of charged particles. 
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When the scattering is independent of the energy and the charge of the particles 
i.e. when R and defined in eq. (|29|) are independent of A and when = S_a, the 
integration over A can be trivially performed as it expresses the completeness of the (fx 
modes. In this case, as in the DF model, one has = [3'^^ = 0. One also finds that 
the emitted flux is simply 

{Tvv) = {Tvv)df , (44) 

where {Tyy) j^p is the flux found in the DF model, see (|TB|). 

Instead, when R and depend on the energy and/of the charge, (3^1^ and (3'^^ will be, 
in general^, different from zero. In this case, one also looses the possibility of re-expressing 
the flux in terms of the derivatives of the trajectory as we did it in eq. (|19|). This can 
be understood from eqs. (^) : when expressing Ta as a series in powers of A, one would 
obtain for {Tvv) a non-local expression in V unless the series in A stops after a finite 
number of terms. 



2.3 Additional remarks 



In this subsection, we relate the matrices Saa' and 8^^^ which act linearly on in and out 
operators to the conventional 5* matrix acting on multi-particle states in the Fock space. 
With this identification we shall be able to relate the Bogoliubov coefficients eq. ( ^Tj ) to 
transition amplitudes and not only to expectation values as in eq. (^3]). 

By definition |24||, the action of this operator on states and operators is the following 



|0m.) — S |0o„t). 



i,out 



s 



s 



I i,out t 



S b 



i,in f 



s 



(45) 



Since we are dealing with elastic scattering, this operator contains exactly the same infor- 
mation as the matrices Sa, Sa. Indeed, the block diagonal character of eq. ( pOD and the 
linearity of eq. (^) tell us that S is the exponential of a quadratic form of proper-time 
operators a\,h\: 



S 



-i[a 



i,in I] 



A 



'AA'^-A' 



7',mt jiAni—ij ij,in\ 



(46) 



Then straightforward algebra gives 



Saa' 
saa' 



6{\ - A') 
saa'(-R, d-, 



Arcsin(i?)e 



Arcsin(/2)e* 



(47) 



where i?, 9, cj) have been defined in eq. (^) and i?, 9, are defined in the same way from 
Sa. We note that in the DF model, i.e. in the limit of perfect reflection, Saa' is given by 



'AA' 



5{\ 




^AA' ' 



(4J 



Although the configurations on the left and and on the right of the mirror completely 
decouple, the S-matrix S treats both sides simultaneously. 



^I3^u ^ requires that T^^ * 



jtC/c/ a > 0, and similarly for the VV coefficients. In the next 

Section, we shall see that condition is satisfied for the time independent couplings with U(l) symmetry. 
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To anticipate the expression of S in terms of Minkowski operators which will mix 
creation and destruction operators, it is a convenient to rewrite eq. (|46|) in term of the 
4- vector a^'™: 



S = e-H«A s^A'«A' ) with {sZ>) = ^J 1 • (49) 



To obtain the expression of S in terms of the Minkowski operators a^'™, 6^'™, it suffices 
to use the matrix B^^x to replace proper-time operators by Minkowski ones. Explicitly one 
obtains 

S = e-H«. ) with sZ> = BZ s^xx' {B^^'x'T" ■ (50) 

Formally, 5* provides the answer to all questions concerning asymptotic states and 
expectation values. For instance, the probability amplitude governing the (Minkowski) 



vacuum decay, eq. ( p.09| ), is simply 

= {Oout\0^n) = {0^n\S\Oin) . (51) 

Similarly, the probability amplitude for an initial quantum of momentum k to be scattered 
and for no pair to be created is 

(0„.,|ar*«r^|0.„) = (0.„|ar S art|o,„) = i(«-i)t . (52) 

The last equality is easily obtained by using eq. (|32D to express a^°"* in terms of a-^'*" and 
^i,o«it. In the same way, the Bogoliubov coefficient /3 is related to probability amplitude 
to find a pair of out quanta in the in vacuum by 

It should be stressed that these relations determine the physical interpretation of the 
overlaps a,P given in eq. (^T]). In fact the second quantized framework was never used 
to obtain eq. (|4l|) : only the linearity of the relations and the orthonormal character of 
the proper time and the Minkowski modes basis were exploited. 

The physical interpretation of a, j3 is the following: to first order in the transfer matrix 
Ta, a (P) divided by Z gives the probability amplitude to scatter a quantum (to produce 
a pair of quanta), since ^ 1 + zT (/5 ~ —iT). Upon considering higher order terms 
in Ta, one looses the simplicity of the relationship so as to get the above equations. The 
simple relation in the linear regime will be nicely confirmed in the next Section, when 
using perturbation theory. We shall see in particular that the division by Z corresponds to 
the usual restriction of keeping only the connected graphs engendered by the development 
of S* = Te~*^^'^*^ in powers of g. We shall further comment on these aspects at the end 
of Section 3. 



3 The self-interacting model 



In this section we introduce a model based on self-interactions which derives from an 
action principle |2^, |10|- In a first part we consider time independent couplings. In 
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this case, re-summing the Bom series leads to diagonal matrices in the proper-time energy 
A with parameters R and (p which depend on A according to the number of derivatives in 
the interaction hamiltonian. This model will be generalized by considering a thick mirror 
with a non-zero width. Using a perturbative approach, we shall see that the thickness 
acts as a UV cutoff. 

In the second part, we work with time dependent couplings. We shall work pertur- 
batively, up to the second order in the interactions. The novelty concerns the transients 
induced by the switching on and off of the coupling. 



3.1 The scattering with g constant 

To exploit the fact that the coupling is r independent, it is convenient to work with the 
coordinates (t, p) in which the mirror is at rest. In these coordinates, the interaction 
Lagrangian reads : 

Unt = g jr^dTdpf{p) J(<l>(r,p),<l>t(r,p)). (54) 



g is the coupling parameter, / is a real function which specifies the thickness of the mirror 
and which is normalized as follows J^^dp f{p) = 1. J is an hermitian operator which 
is quadratic in the complex field. We shall consider three different cases: + 

^Hd-r^, and dr^^dr^ + dr^dr^^ ■ In the following equations, we shall present the details 
only with the second expression. At the end of the derivation, we shall give the final 
results for the two other cases. 

Given eq. (|5lD, eq. (|I]) is now replaced by 

id'^-d'^Mr,p)=gfip)2tdr^. (55) 

Being linear, the solution can be expressed as 

rr+oo 

^T, p) = <|.-(r, p)+g dr'dp' G^-^\t, p ; r', p') /(p') 2t ^,,^T', p') 



oo 



<|.-*(r, p) + g dr'dp' G^'^^(r, p ; r', p') /(p') 2i 9.,<|.(r', p') , (56) 



oo 



in terms of the homogeneous solution ( <|>°"*) which determines the initial (final) data. 
The retarded and advanced Green functions are defined, as usual, by 

G"'(r, p ; r'. p') = //_ .A ,1 - _ „ ^ ,,,, (= fa r' > r) , 

G-(r. p ; r', p') = //_^ „ (= for r' < r) . (57) 

To exploit the time independence of the coupling g, we work at fixed energy with 

^a(p)=/ dr — ^T,p)e''-. (58) 

J-oo ZTT 

In Fourier transform, eqs. (|56|) give 

mp) = ^Tip)+^9 r^dp' f{p')^,{p')e^\''-'''\ 

J — oo 

P + OO 

= <ip)-W dp' f{p')v,{p')e-^'\^-P\ . (59) 
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These equations have been obtained by using 



hoo 



P^'-(P-P') +99V 



00 P-{\±ieY 2{\±ie)'' • ^^^^ 

We now decompose the quantized modes in terms of creation and destruction 
operators, 

(^r(p) = ^^{a^'re^P + aY-e-^'^^) for A > 

= ( 6[lf ^ e-^l^l" + fef^^f ^ e'l^l^ ) for A < 0. (61) 
47r|A| 

We do the same for the out modes. Then, for /(p) = 5{p — po), in the hmit e ^ 0, eqs. 
(H) give : 



(62) 



We recover the hnear structure of Sa in eq. (p5|). Since the unitarity of provides 
S'a * a^'*", when using the definitions of eq. (EW), we obtain 



R = I , = Arctan(5f) and Q = 2Apo • (63) 

vi + r 



In the strong couphng hmit {i.e. for g +00) one obtains total reflexion (PB|) in a A 



independent manner. This is a special feature of the coupling J = $^0t-$ which is 
associated with a dimensionless g. 

This analysis can be repeated with the two other operators previously defined. The 
presence or the absence of derivatives in J modifies the IR or UV behavior of R. For 
$t$ _^ <|)$t one obtains : 

Rx = . , 0A = Arctan(^/A) and 9 = 2Apo • (64) 

'1 + gyx^ 



In this case, the mirror is totally reflecting in the IR. This leads to strong IR divergences 
when considering time dependent coupling g. On the contrary, when using dr^^dr^ + 
dr^dr^\ we get 

-^^ = 5 (f>x = Arctan(5(A) and 6 = 2Apo • (65) 

Vl + g ^ 

In this case, the mirror is transparent in the IR limit. This useful property will be 
exploited in Section 4. 

One notices that the transfer matrix can be expressed in a general way according to 
the number n of derivatives d-r in the interaction term: 



1 - igX'^-^A, [ 1 



(66) 
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In this expression, we have not taken the hmit e ^ 0. The function = A/(A + ie) 
determines the analytical properties of Ta in the complex A plane. The specification of 
the pole of follows from that of G^'^* in eq. ([571) . It guarantees that causality will be 
respected p6|. This crucial ingredient was missing in Section 2.2 wherein the matrix Ta 



can be chosen from the outset. In that kinematic framework, the analytical properties 
should be imposed by hand if one wishes to implement causality. On the contrary, in the 
present case causality follows from the Heisenberg equations. 



Eqs. (^) and (61) also determine the relation between the anti-particle in and out 
operators b\ . By direct computation one finds = — T„a- This is precisely the condition 
which gives = P'^^ = 0, see the former footnote. When using Ta, Ta in eqs. (|41|), we 
obtain the Bogoliubov coefficients relating inertial modes. And from these coefficients. 



one gets the mean value of energy flux Tyv as in eq. (|^), but with causality built in. 

We now study the case of a thick mirror with J given by ^Hdr^. To display the 
effects of /(p), it is convenient to work with the (spatial) Fourier components. Eqs. ([59| ) 
become 



2g\ 

(A + ze)2 -/2 



+ CXD 



OO 
+ 00 



For an arbitrary window function /, these equations do not lead to analytic relations 
between asymptotic in and out fields. Therefore, we use perturbation theory. To first 
order in g we get 




For a normalized gaussian function / centered on Pq, the non-diagonal terms which de- 
termine the reflection probability are ge^'^^^P'^e~'^^ ^ . Therefore, a, the spread of the 
mirror reduces the reflection of high frequencies: for X ^ 1/a the mirror is completely 
transparent (this is also true for the two other J's). 



3.2 The scattering with g time dependent 

In this subsection, the coupling parameter is a function of the proper time (?(r) = g f{T) 
where /(r) is normalized by J^^dr fij) = 2T with 2T the proper-time lapse during 
which the interactions are turned on. Unlike what we had in the former subsection, 
resumming the Born series is no longer possible since the time dependence of the coupling 
destroys the decoupling of the equations into sectors at flxed frequency A. In fact we meet 
a situation analogous of that of a thick mirror which mixed different momenta. Thus, we 
shall work perturbatively: all quantities will be evaluated up to the second order in g. 

We remind the reader that in the interacting picture, the operator $ evolves freely, 
i.e. with g = 0. Therefore the in operators a^, 5^ coincide with the out operators and are 
equal to the usual Minkowski operators. Hence, they deflne the (Minkowski) vacuum |0). 
Instead, the states evolve through the action of the time ordered operator: 

|^(t = +oo)) = Te^^ |^(t = -oo)), (69) 
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where L = g J dTf{T)J. Since the trajectory is timehke the time ordering with respect to 
the Minkowski time t is equivalent to that of the proper-time r. 

To make contact with Section 2, we shall work in this section with the state \^o{t)) 
which is equal to |0) for t = r = — oo. When expressing its final value in the basis of the 
unperturbed states, i.e. the states which would have been stationary in the absence of 
interactions, we get 

\Mr = +oo)) = |0) + E// dcoduj'{BZ,+CZ>)\uJuj')^^ (70) 



where 



^1' = W (0| al li, [j^yr f{r) J(r)) |0) , \ujuj')^^ ^ a}J lPj\Q) , (71) 

dr / dr' f{r) f{r') .J{r) .J{r')) |0),. (72) 



We have limited the expansion in g to these three terms since we shall compute the 
energy-momentum tensor up to g^ terms only. As before, i,j denote the U,V sectors 
and uj,uj' Minkowski energies. The symbol ( )c means that only the connected part of 
the expectation value is kept. This restriction follows from the fact that the contribution 
of the disconnected graphs cancels out since they also appear in the denominator of the 
expectation values, see e.g. ||24 |. 

Using eq. (|70D, the expectation value of Tyv is given by : 



(Tyv) = (^o(r = +oo)|Tvy|^o(r = +oo)) 



3?e 



- e-^(-'^-) (i^rj + i^rj + c^rj + , (73) 



where B^ji^ and C^J, are related to the unbar quantities by inverting particle and anti- 
particle operators, thus B'J^, = B^]^ and CjJ = C^,^. 

Since the integral of the second term in eq. (73) vanishes and since bar and unbar 
quantities differ at most by a phase, the total energy received on the V part of J'^ is : 

/+00 /"OO /"OO 

dV{Tvv), = 2Y. duu dk\Bll\\ (74) 
-oo j JO JO 

Hence only the B terms contribute to the energy as the P terms did in eq. (|T7|). 

In order to compute the local properties of the flux, we need to compute the second 
term of eq. (|T3|). To this end we decompose C^^, into two parts : 

c^rj=i?L^-(oiarfer'^^io),, (75) 

where 

R:^' = 'Imalbl, LL |0),, (76) 

V = ^(Prfr r^dr' f{r)f{r')e{r-r')J{r)J{r')) (77) 
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and e(r — r') = 9{t — r') — 9{t' — r). Then, {Tvv)ti, the contribution of V to (Tyy), 
enjoys the following properties. First it carries no energy. This is obvious since it is 
built with terms which all contain e*^('^+'^). Secondly, it vanishes for /(r) = est. This 
can be understood from the fact that the time ordering properties can be encoded in the 
analytical properties of the matrix which is diagonal in A, see eq. This means 

that this term modifies the shape of the transients related to the switching on and off of 
the interaction but without affecting their energy content. In the rest of the paper, we 
shall therefore ignore this term. 

We now compute R^^>- Since only the connected part should be kept, we can insert 
the following operator between the two operators L in eq. (|76D 

Y: Hdk' al^ bi^\0){0\ K' ■ (78) 

Grouping together, as in eq. ([73|), the first order and the second order contribution in g, 
we get 

_ _ pOG _ _ 

+ + RZ^> + i?L^ = T. dk {A^l * B:,i + B^l) (79) 

■ JO 

with 

= (0| al (l + zL) ai^ |0), and = (0| bl (l + zL) bi^ |0), . (80) 

Hence we find that {Tyy) is given by eq. (E2|) with 



^r^^^yart^cle ^ T du H du' (81) 

10 Jo 2n 
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[e-^(-'-) {[dkB:i:B:?,)-e-^^-'--^ {fdkAirB:^) 

(Tvv) anti-particle gi^en by the same expression with All,B^l replaced by All,B^l. 

Thus, to second order in g, we recover the structure of eq. (^) which gives the flux 
emitted by a partially transmitting mirror. The Bogoliubov coefficients a^l and P^l 
have been replaced by the transition amplitudes A^l and B^^. In this we recover the 
correspondence of eqs. ( |52[ ) and (^) when considered to first order in the transfer matrix 
T^. This is not surprising since the evolution operator Te*^ which defines A*^ and B^^j^, 
given in eqs. ( |80| ) and ([71]), is the operator S of eq. (|i5|). 

This correspondence is nicely illustrated in the case where g{T) = g and J = ^Hdr^- 
In this case, to order g but whatever is the mirror's trajectory U = Ud{V), one has the 
following identities 

Al^ = ga^k, B^t' = 9f3.k, (82) 

where a^k and are the Bogoliubov computed in the DF model. These relations 
establish that a^k and P^^k should be understood as transition amplitudes. This implies 
in particular that the momentum transfers to mirror (which have been neglected so far) are 
respectively h{k+uj) and h{~k+uj). This imposes limitations when considering ultra-high 
(trans-Planckian) frequencies since neglecting the momentum transfers requires hu -C M 



where M is the mass of the mirror |T0[. Thus, when high frequency quanta are emitted. 



the validity of the predictions obtained with a recoil-less model must be questioned [0]. 
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4 Applications 



First, we analyze the properties of the transients associated with the switching on and 
off when the mirror is at rest {z = 0) and in Minkowski vacuum. Secondly, we generalize 
this analysis by replacing the Minkowski vacuum by a thermal bath. Then, we use the 
well-known parallel between inertial systems in a thermal bath and uniformly accelerated 
systems in vacuum to obtain a regularized expression of the flux emitted by a uniformly 
accelerated mirror. 



4.1 The transients in vacuum 

We first focus on the frequency content of the transients. For an inertial mirror at rest 
at 2; = in Minkowski vacuum, the transition amplitudes A and B of eqs. ( |80D and ( [7I| ) 
can be expressed in terms of the Fourier transforms of f{t) 

U = ^ jdt f{t)e'-' . (83) 



To order g we obtain 



= Siu-k)S'^+zgU_,^-^ , (84) 



'Lok 



where 

{1 for <l>^$ + 

u + k for <^Udt<^ ■ (86) 

uk for dt^^dt^ + dt^dt^^ 

Thus to order g'^, the mean number of V particles of energy uj is given by : 

W) = 2 E / \B^,\' . (87) 

The factor of 2 arises from the fact that the pair production amplitudes B are independent 
of the (relative) sign of the momentum of the two quanta. 

Our aim is to describe the transients associated with the switching on and off of the 
coupling to the mirror. To this end and to be concrete, we shall work with the function 

m = i(tanh(^) - tanh(^)) . (88) 

It is constant during a lapse of time 2T centered about t = and the time intervals of 
the switching on and off are ~ 4 A. In the limit A — 0, / tends to the square window 
[6{t + T) — 6(t — T)]/2. The Fourier components of / are 

f = - ^^^^^^^ (89) 
2 sinh(a;7rA/2) ' ^ ^ 
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One sees that the UV behavior is exponentially damped by A. On the contrary, in the 
IR, — > T/tt since the coupling lasts during 2T. 

When considering the first two cases of j{uj, k) of eq. (|86|), these simple observations 
imply that the mean number (N^^) is ill defined since the integral over k in eq. (|87| ) 



diverges in the IR. Therefore to obtain well defined expressions we shall work with the 
third case. In this case one has 

,,,,,, g^A^ r,, , iAn^{{uj + uj')T) , , 

W) = ^^/ d^'^' ■ ,2,y ^ J.L • (90) 
/ JO smh ((w + tu'jTrA/i) 

It is perhaps appropriate to discuss the condition on the (dimensionfull) coupling 
constant g which guarantees the validity of a perturbative treatment limited to order g^ . 
The condition is that the mean number of quanta per quantum cell (which is equal to 
{Nj)duj ^ {N^)n/T in the limit uT ^ 1) be well approximated by eq. (|87D . This requires 
that the probability to obtain two quanta in a cell is negligible with respect to that to 
obtain one. This translates mathematically hj g'^ ^ T A in the limit of interest T/A ^ 1, 
i.e. when the flat plateau is much longer that the slopes. The condition g^ <^T A means 
that the limit T — > cx) can be safely taken. Instead the limit A — is more delicate. A 
sufficient condition consists in working at fixed g^ <ti 1 where g = g{T A)^^^"^ . A stronger 
condition is to impose that the total number of particle emitted, /q°° {N^)du!, is finite in 
the limit A 0. In this case, g = g/A should be held fixed. 

When studying eq. (pO]), one first notices that in the limit T ^ oo with g and A fixed 
the total number of particles emitted is independent of T, thereby not giving rise to a 
Golden Rule behavior characterized by a linear growth in T. Secondly, (A^^) is maximum 
for oc 1/A. Finally, for uA ^ 1 one has (N^) — e~'^^^. We thus find all the expected 
attributes of transients: their particle content is independent of the duration T when 
T/A ^ 1, and their Fourier content is peaked around the adiabatic switching rate A~^. 

We now study the space time repartition of the energy fluxes associated with these 
transients effects. We first notice that once the P-term defined in eq. (^) as been 
dropped, the mean flux can be expressed as : 

(Tyy) = -253m ((OlTyyLlO)) + 3?e ((0|L [Tyy, L] |0)) . (91) 

Of course, by decomposing L and Tyy in terms of creation and annihilation operators, one 
would recover respectively the linear and the quadratic contribution of eq. (0). However, 
being interested in the space time properties, we shall not do so and shall work instead in 
the time 'representation' with the V part of the (positive frequency) Wightman function. 
It obeys 

dvWiV - V) = dvmKV, U) ^V, U'm = ] . . (92) 

Using this function, the linear contribution in g reads 



dtf{t) {d^dvW{v-t)} 



127r 



dU{t = V) . (93) 



To obtain this result, we have integrated by part three times. The boundary contributions 
all vanish since / given in eq. (^8]) decreases faster that any power of t. The last integration 
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is trivially performed by using 53m((x — ie)^^) = tt 6{x). These properties explain the 
local character of the expectation value^ 

To evaluate the quadratic contribution in g we proceed along the same lines. We first 
evaluate the commutator so as to obtain a quadratic form in $ and $y, where $y means 
only the V part of the field operator $ should be kept. We notice that the derivatives dt 
in J might be expressed as dy since they are evaluated at z = but they act both on the 
V and the U part of $. Using this notation, one finds 



[Tvv,L] 



ig fiV) [{dv^l 4$ + dv^v 5^$^) + (h-c 

(h.c 



Then the g"^ contribution of Tyv is 

{Tvv),uadr = fiy) (t J dV f{V') {dy,d'yW{V' - V)) {dv'dyWiV - V)) 

+16/ dvf (^i J dV f{y') {dv'dvWiV' - V)) 



(94) 



12tx 



{fd'yf + 2dvfdlf 



(95) 



Since neither / appears in eq. ( p3| ) nor p in eq. (0), one recovers the fact that an 
inertial mirror doesn't radiate while its coupling is constant. This is illustrated in Fig. 1. 
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Figure 1: The thin soUd Hne is J{t) given by eq. for T = 10 and A = 1. The dashed 

line is {Tvv) quadr ^'^'^ ^^"^ thick line is the part of (Tyv) which contributes to the energy, see eq. 
(^. These two curves have been plotted in the same arbitrary units. The behavior of {Tvv)iin 
is similar to {Tvv)quadr- 

^ It should be pointed out that we could have written {Tvv)iin b& a commutator. This however is 
not appropriate since one looses the analytical properties of W which are encoded by ie (They arise from 
frequency content of the vacuum and play a crucial role in defining the above expressions) . By performing 
first the commutator (or equivalently by first taking the imaginary part in eq. (|9^)) one would obtain 
an ill defined expression. The same remark applies to the quadratic term in g. To obtain well defined 
expressions, only one commutation (and not two) should be done. 
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As in eq. ([T9|), one can decompose {Tyy) into two parts : a positive definite term and 
a total derivative which does not contribute to the total energy 



(Tvv) - {Tvv)iin + {Tvv) 



9' 



Thus the total energy is 
(Hv) 



12tx 



quadr 



(96) 



dV (dlf) 



127C J-oo \ / JO 



dujuiN::;). 



{Hv) is finite when the mean number (A^o;) decreases faster than uo ^. 
when working with eq. (|88D at fixed A 7^ 0. In this case, one finds 



= ir. AS a) ■ 



(97) 

This is the case 
(98) 



The main feature {Hy) is that it is independent of T in the limit T/A ^ 1, see Fig. 2. 




Figure 2: The plot Fi{x) defined in eq. 
X = T/A. 



|) in such a way that F ^ 1 for a; — > 00 where 



m 



In conclusion, we consider the limit A ^ which corresponds to the situation studied 
in view of its analogies with the residual flux emitted at the end of the 
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and in 

evaporation of a black hole. In this limit, f{t) becomes a step function, the energy flux 
is concentrated in a narrow lapse A and its frequency content diverges. In fact (Tyy) 
becomes a distribution since it is built on the derivatives of f(t). However the singularity 
is worse than a delta, as clearly seen from eq. (^5)) . This means that the total energy 
emitted is also singular, as indicated in eq. (PB]). Moreover, there is no reason for the 
nature of the singularity to be universal. It depends on the number of derivatives in 
the Hamiltonian and it might also vary when considering higher orders in g. Hence the 
question: what is the flux emitted by the disappearance of the reflection condition, is 
not univocally defined. To have a well-defined question, one should first choose a regular 
model such as that defined by eq. ( p^ with J containing at least two derivatives and g{t) 
given by eq. (|88|), and only then consider the singular limit. 
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4.2 Flux and energy in a thermal bath 

Our aim is to obtain a regularized expression for the flux emitted by a uniformly accel- 
erated mirror. To this end, we shall use the isomorphism between the flux emitted by 
a mirror at rest in a heat bath at temperature a/27r and the Rindler flux emitted by a 
uniformly accelerated mirror of acceleration a. 

In a thermal bath, the V part of the Wightman function obeys 



dvW^iV - V) = --^^ coth { ^{V - V 
Ait p \p 



. (99) 

It reduces to dyW of eq. (|92[) in the zero temperature limit, i.e. for (3 oo. When 
replacing W by IV^ in eqs. (^) and (^) we obtain the mean flux emitted in a thermal 
bath. It can be shown to beQ 

{Tvvf = ^dlf-^(f{y)d'yf + 2dyfdlf 

i/TT iZTT ^ 



9 ^ ^ 9^ 



^^^dvf - ^ (fiV) dlf + 2dyf dvf) ] . (100) 



The first two terms are equal to eq. (|96[) and the last two scale like (A//?)^. Thus, they 
are negligible in a low temperature limit, /? ^ A, and dominant in the high temperature 
regime. 

We are now in position to obtain a regular expression for the flux emitted by a uni- 
formly accelerated mirror in Minkowski vacuum. Using the well-known isomorphism be- 
tween systems at rest in a thermal bath and accelerated systems in vacuum, the mean 
flux of Rindler energy emitted by a mirror of acceleration a is 

{TMr^ = {TvviV = v)f , (101) 

where v is the null advanced Rindler time {av = In(al^)) when the mirror is located in 
the Right Rindler quadrant, {z < \t\). The coupling between the mirror and the field is 
turned on during a proper time lapse 2T and the switching on and off rate A~^ is also 
measured with the proper time. 

In the limit T ^ (3 and a~^, (T„„(f ))°'^^ ^ at fixed |f | < T since the flux is localized 
in the transients of 'thickness' A. In this we recover the well known result that a uniformly 
accelerated mirror does not radiate. In the DF model, this immediately follows from eq. 
(P^). (Notice that this vanishing is a universal property of accelerated systems when they 
have reached equilibrium with the Rindler bath|2^, ^ |1^.) However this vanishing 



flux is accompanied by transients effects whose Minkowski properties become singular in 
the limit T —>■ oo whatever is the value of A. To have regular Minkowski properties means 
that the total energy 

/ + 00 
dv {TZ) , (102) 

and the mean number of Minkowski quanta emitted by the mirror be finite. 

When requiring that {Hy'^) is finite, the coupling /(r) must decrease faster than the 
Doppler effect dV/dv = e""^ . Using eq. (^), this implies a A < 2. (A similar condition 

^The details of the calculation will be presented in p3]. 
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also applies when considering the fluxes emitted by an accelerated two level atom|]T6|. 
Using eqs. ([lOOD , ([lOlD and ([102D , we get 



127r 



dv e° 



(103) 



What we have learned from the condition aA < 2 is that the Rindler fluxes which lead 
to finite Minkowski energies are dominated by the vacuum effects governed by A rather 
than by the temperature effects induced by the acceleration. 



A The in — out overlap in the non-stationary case 

In order to have simple expressions for this overlap, we will use a discretized basis of wave 
packets in which the integrals are replaced by sums and Dirac distributions by Kronecker 
symbols. 

Instead of working with the in and out basis, it is appropriate to define a third class 
of operators d^,b^. This new basis generalizes the Unruh modes ||25|, |[ since (6^) is 
made out of a^" (6^") but is characterized by a fixed out frequency u: 



0^f,i CLi, 



Y.<kaT, aJ^ = Y.al,U^. (104) 



The real coefficients a^^ are such that [d^,dj^] = 1, therefore = Sfc |«a;A:P- The notion 
of particle/anti-particle is obviously maintained since the a are made of a*" only. Notwith- 
standing, for arbitrary a^^k and fi^^k this new basis is not orthogonal and the commutation 
rules are given by 

[d., dl] = F^^, = [i, bl] = ^ . (105) 

By construction from eq. (0), these new operators are related the out operators by 



out_^ I u With 5,,, ^ E/^-'^^L'', (106) 



where a^J is the inverse matrix of a^k- ( «wfc is always invertible since otherwise there 
would exist incoming particles whose scattering would give only anti-particles.) 

As for the Unruh modes, the new basis is useful to relate in a simple way the out 
vacuum to the in vacuum. Straightforward algebra indeed gives 

Kut) = I exp ( E ^ F-k Bk.' dill) |0,„) , (107) 

where Z is defined by 

= \{Oout\0^n)\^ . (108) 

Even though eq. ( |107| ) looks cumbersome, one easily verifies that, to order it 
correctly gives the relationship between the vacuum decay {Z > 1) and the pair creation 
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probability of Minkowski quanta. Indeed, using {F ^)ujw'{0in\O'u}' o^'/|0i„) = S^^uj" and the 
condition on B^^/ and F^^/ which arises from [a""*, b^^] = and eq. (|105|) , one obtains, 

Z2 = l + ^ \b^^,\\0{P^) . (109) 

UlUj' 

This is the correct expression since the probability to have a pair of out quanta is 

|(Oo„t|ar &°^*|0.„)P = = \B^^'\' + 0{P') (110) 

For completeness, we notice that when the scattering is stationary (as it is the case 
for uniform acceleration and in black hole evaporation), one has 



Since they are diagonal, eq. (|107|) becomes 



\0out) = ^exp{J2 — ^lbl) \0^n) (112) 



thereby recovering the usual diagonal expression governed by the "Unruh" operators 
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